Quantum critical scaling in graphene 
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We show that the emergent relativistic symmetry of electrons in graphene near its quantum critical 
point (QCP) implies a crucial importance of the Coulomb interaction. We derive scaling laws, valid 
near the QCP, that dictate the nontrivial magnetic and charge response of interacting graphene. 
Our analysis yields numerous predictions for how the Coulomb interaction will be manifested in 
experimental observables such as the diamagnetic response and electronic compressibility. 



Recent experimental developments have made possi- 
ble the study of graphene, a single-atom thick sheet of 
graphite [lj]. The novel electronic properties of graphene 
arise from the linear, cone-shaped energy-momentum dis- 
persion of electrons at low energies. This condensed mat- 
ter realization of a relativistic Dirac spectrum follows 
from simple models of electrons hopping on the honey- 
comb lattice of graphene|2L and has been confirmed by 
a range of experiments [EUl Q ■ 

The relevant Hamiltonian is that of relativistic 
Coulomb-interacting fermions in two dimensions: 
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with velocity v ~ 10 8 cm/s[3|. p; = ~ihW ri is the mo- 
mentum operator and a — (a x ,a y ) are Pauli matrices 
that act in the space of the two sub-lattices of the honey- 
comb lattice structure. There is an additional N = 4 fold 
degeneracy caused by spin and the two distinct nodes of 
the dispersion, with the only tunable parameter in Eq. ([1]) 
being the dielectric constant e. 

In case of the usual electron gas, the first term in 
Eq. (JTJ) is J2iPf I '(2m). Then, dimensional arguments 
imply that the kinetic energy dominates for high elec- 
tron density while the Coulomb interaction dominates at 
low density. The linear Dirac spectrum changes this sit- 
uation. The relative importance of the potential and ki- 
netic energy is the same for all densities and controlled by 
the dimensionless number A = e 2 / '{AevTi) . For Acl, the 
Coulomb interaction is negligible. Using the above value 
for the electron velocity yields A ~ 0.55/e, i.e. A ~ 0.55 
for a free standing graphene film in vacuum, implying 
that one cannot ignore the Coulomb interaction. The 
role of interactions in g raphene has been discussed pre- 
viously [1, 0, 0, 0, Ea El El EE El- However, few 
specific predictions for observable quantities have been 
made that allow for a comparison with experiment (see, 
however, Ref. [Hj]). Here, we exploit the enlarged sym- 
metry near its quantum critical point (QCP) to deduce 
numerous predictions (based on scaling theory) of inter- 
acting graphene. 

In this letter, we use a renormalization group (RG) 
approach to the Hamiltonian, Eq. |T]), and analyze the 
magnetic and charge response of graphene as a function 
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(1) FIG. 1: (Color Online) Quantum critical phase diagram of 
graphene as a function of density n (in units of 10 12 m -2 ) and 
temperature T (in K), for the vacuum case e = 1, showing 
the Dirac liquid and Fermi liquid regimes separated by the 
crossover temperature T* [dashed lines, Eq. (I13|l ]. with the 
quantum critical point occuring at n — T — 0. 

of temperature T, carrier density n, chemical potential /i 
and magnetic field B. We make specific predictions for 
the the compressibility k — dn/dfj,, the diamagnetic sus- 
ceptibility xd, the magnetic moment M(B), the heat ca- 
pacity C, the infrared conductivity u(uj) and the density- 
density correlation functions x c (q,u;). We demonstrate 
that interaction effects in these quantities are measur- 
able, allowing experiments to reveal the subtle interplay 
of interactions and kinetic energy in a Dirac liquid. Our 
analysis is based on the fact that for T = B = /j, = n = 0, 
clean graphene is located at a QCP, as illustrated in 
Fig. H and its properties nearby can be obtained via 
crossover scaling arguments. 

The low energy action that follows from Eq. (TT]) is 



S = % / ^ (d T ao — M»V r -a) ip + — 



(2) 



Here, ijj = tp (x) is a two component electron field where 
x = (r,T, s) stands for the 2D position r, imaginary time 
r, and valley and spin quantum numbers s — 1 • • • 4, such 

that J x ... = J d 2 r f£ drY,t=i • ■ ■ with Z 3 " 1 = k B T/h. 
n-x — (%) ^ (x) is the electron density. 

We perform a one-loop Kadanoff- Wilson RG analysis 
of Eq. ([2|). Fourier transforming ip(x) yields tp(k) where 
k = (k,w n , s) with planar wave vector k and Matsubara 
frequency u n = (2n+l)kBT/h. We trace out high energy 
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modes with A/6 < |k| < A and obtain a renormalized 
action. Due to the static nature of the Coulomb interac- 
tion, the one loop fermion self energy X (k) is frequency 
independent, i.e. the fermion dynamics remains unrenor- 
malized: d T — * d T . Higher order diagrams @, [l(| or strong 
coupling effects^ cause X (fc) to be w-dependent but are 
beyond the leading diverging terms considered here. Sim- 
ilarly, vertex corrections describing interactions between 
electrons and collective charge fluctuations vanish. This 
is consistent with the Ward identity that follows from the 
conservation of the total charge and implies e 2 — > e 2 un- 
der renormalization. The only nontrivial renormalization 
is that of the velocity. Since 9E (k, to = 0) /9k diverges 
logarithmically as |k| — > 0, within the RG approach this 
yields v — > u(l + Alog6). To complete the RG proce- 
dure, we rescale the fermion field 



V>(p,w,s) -> (b)ip(bp, Zrj, 1 (b)u,s), 



(3) 



with Z T (6) = &- 1 (l + Alog&) and (6) = frZ^ 1 (6) 
that define the relationships between fields in the original 
and renormalized theory. The increase of the velocity 
yields a decrease of the dimensionless coupling A. 

Upon iterating the RG transformation one finds the 
RG equations for A (b) as well as the temperature T (b): 



d\(b) 
din b 
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dT(b) 
dln& 



T(6)(l-A(6)), (4) 



These equations are solved by T(b) = Z^ 1 (b) T and 
A (6) = Xb~ 1 Z^ 1 (&). The coupling constant is marginally 
irrelevant and will lead to logarithmic corrections relative 
to the Dirac gas of non-interacting electrons with linear 
spectrum. Just like in other quantum critical phenom- 
ena, the temperature is a relevant perturbation, causing 
crossover to a classical, finite T regime. Interactions lead 
to T (b) < bT, i.e. the quantum dynamics of the interact- 
ing Dirac liquid is more robust against thermal fluctua- 
tions than the non-interacting Dirac gas. The coupling- 
constant flow equation was obtained earlier in Ref. [6[. 

We now utilize the RG to develop general scaling rela- 
tions for physical observables. We start with the electron 
density n and the compressibility n = dn/dfi. A finite 
compressibility can be the result of a finite chemical po- 
tential fj, or due to thermal excitations at fj, = 0. To 
analyze the system for finite chemical potential we add a 
term — fj, J rp'aoip to the action and rescale the fermion 
fields yielding fx(b) — Z^ 1 (b) ii. Upon renormalization, 
the number of particles per area obeys: 

n (T, fi, A) = b- 2 n {Z T X (b) T, Z T X (b) fx, X (b)) . (5) 

This equation may be written in the short-hand nota- 
tion n = b~ 2 n,R, which we shall use henceforth, with 
the subscript R denoting renormalized quantities. Thus, 
performing the derivative with respect to /i yields for the 
scaling of the compressibility k = b~ 2 Z^ 1 kr. As nu- 
merous similar results appear below, we now discuss the 
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FIG. 2: (Color Online) (a) Plot of the inverse compressibil- 
ity k -1 = (in units of (10~ 10 meVcm 2 )) as a function of 
density n (in units of 10 12 /cm 2 for the case of v = W 6 m/s 
with the dielectric constant e = 5.5 and e = 1 (solid) and 
the noninteracting case (dashed), along with the data from 
Ref. (b) Plot of the magnetization per area (in units of 

fj,A) as a function of magnetic field (T) . 

physical meaning of this expression in some detail. The 
left side of Eq. (JSJ) is the physical density in graphene, 
which is the quantity of interest. The right side is the 
density in the renormalized system where the effective 
coupling A(6) is small and the effective temperature T(b) 
is high. With an appropriate, physically-motivated choice 
for the renormalization scale b, we can put the renor- 
malized theory into a regime in which the calculation is 
particularly simple. We fix b by noting that the RG 
equations were derived assuming the low-T (quantum) 
limit, and are thus only valid for T(b) < Tq = DJks 
with the bandwidth D — KvA. Thus, we choose the 
renormalization condition T(b*) = Tn [17| which yields 
b* = (l + A log Tq/T. Approximating the renormal- 
ized high temperature compressibility by its free fermion 
result, i.e. kJ, 1 ~ ir (h~v) 2 / (4fcB7oln2) we obtain 
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describing the nontrivial temperature dependence of the 
compressibility of graphene, valid in the shaded region 
of Fig. [TJ Since n (T, fx = Q) = 0, we can determine the 
density as a function of T for finite chemical potential 
/i <C fceT as n(T, A, /i) ~ k (T) \i. At T = but finite /x 
we use n(b*) — hvA and it follows in full analogy 
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and k (fx) = dn/dfi ~ 2n(/x)//i. Since it has re- 
cently been measured in scanning single electron tran- 
sistor experiments flij] we also report the compressibility 
as function of density (with hq = A 2 /ir): 



K_1 {U) = UV \[^\ ( 



1+ 2 log R 



(8) 
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that was also recently obtained by Hwang et al [151 ] . 

For a quantitative analysis we need to make a choice 
for the upper cut off. Following essentially Ref. [HI we 
choose A such that A 2 = 2tt/A where A a = 3 3 / 2 a 2 l /2 
is the area of the hexagonal unit cell with C-C distance 
a = 1.42 x lCT 10 m. It follows T ~ 8.34 x 10 4 K, D ~ 
7.24 eV and n ^ 3.8 x 10 15 cm~ 2 . Thus, the arguments 
of the logarithms are very large for realistic parameters 
of T, /i or n, making correlation effects important. 

In Fig. [2k- we compare experimental data for k^ 1 (n) 
of Ref. [16| with Eq. (jSJl for different dielectric constants. 
For e = 5.5, excellent agreement between theory and ex- 
periment is obtained. While the data can be fit by the 
Dirac gas expression Kq 1 = hv^Jj^ with v as a fitting 

parameter, it cannot be understood without incorporat- 
ing interactions, since the measured velocity v is known 
quite accurately from other experiments. To observe the 
logarithmic variation of k with n, a larger density regime 
or a combined T-n scaling analysis would be required. 

We continue our analysis with the diamagnetic sus- 
ceptibility xd — §jf OI graphene, where B is the mag- 
netic field and M is the magnetization per area. It was 
previously shown [l9|, [2(| that, for a Dirac gas, xd di- 
verges at T — B — 0. We will show that interac- 
tions enhance this divergence. To determine the scal- 
ing properties of \d, we add a gauge-field A fJ- (x), to 
the action: S — ► S + J x A^(x)j^ (x) with the electrical 
current j t _ t (x) — evtfj^ (x)cr^(x). By repeating the RG 
analysis in the presence of A^(x), we obtain the scal- 
ing properties for the Fourier transform of the current 
j M (k, uj) = ZjZ^j^bk, Z^ui) that relates the physical 
current to the current in the renormalized system. The 
factor Zj = X/X(b) follows from a nontrivial vertex cor- 
rection for the current operator and reflects the fact that 
jfj,(x) is a nontrivial composite operator (i.e., it is com- 
posed of more than one if) field) that mixes high-and-low 
momentum degrees of freedom. For the diamagnetic re- 
sponse, we are particularly interested in A„(r) that repre- 
sents a magnetic field (pointing in the z direction, normal 
to the graphene sheet) via B = z ■ (V X A). Standard 
linear-response arguments dictate xd is related to equi- 
librium fluctuations of the current: 
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■K xy (q,0), 



(9) 



where -£f a/3 (q,w) = (j Q (q, w)^(-q, u)) is the current- 
current correlator. The scaling relation K a p — ZxK a p r 
follows directly from that of j Q (k, to) yielding xd — 
b 2 ZtXd,r, and, following the scaling analysis of k: 



^ (T) ^6^( 1 + Al0g Y 



(10) 



where we used the diamagnetic susceptibility of non- 
interacting Dirac fermions[20|. Thus, the RG analysis 



has revealed a nontrivial enhancement of the diamag- 
netic response (relative to the non-interacting result) of 
intrinsic graphene, a result that we have verified pertur- 
batively to order A, following the procedure of Ref. [2l[ . 

Our approach sheds light on the so-called universal 
finite-frequency conductivity 22] (t(uj) of graphene, re- 
lated to the current-current correlator K a p via the Kubo 
formula a(oj) = ^-ImiQ^ (q= 0,w), where the super- 
script ret. indicates the retarded function obtained via 
iu> n — > lu + i0 + . Although Ar Q( g(q, uS) has nontrivial scal- 
ing behavior, the Kubo formula yields a — &r for a(u>). 
This simple scaling relation is constrained by a Ward 
identity to be valid to all orders in perturbation theory, 
as can be shown using an argument due to Gross [231 ] . 
and implies that the T — > 0, w — > (with lu > T) uni- 
versal conductivity is independent of interactions. To see 
this, we take (as we have above) the perturbative re- 
sult for the right side, yielding a{uj) = a(Z^ r 1 cj) = 
Thus, we find that interactions do not modify the uni- 
versal conductivity, in disagreement with recent results of 
Mishchenko Again, this result may be verified per- 
turbatively by computing the leading-order corrections 
to the non-interacting universal conductivity; the crucial 
divergent parts of the leading diagrams identically cancel. 

Our results for the scaling properties of xd, k etc. can 
be alternatively derived from the scaling properties of 
the free energy energy density F(T, A, /i, B). It holds 
F = b~ 2 ZTFn, a result that is easy to verify to leading 
order in perturbation theory; more generally it may be 



derived using the method of Ref. [24(. The scaling of 
the external magnetic field B follows from the scaling 
of the gauge field A. Since within the RG we rescale 
positions via r — * r' = r/b, and A enters via minimal 
coupling, we must have for the latter A — > A' = 6A and 
therefore (from B = V x A) B (b) = b 2 B. This allows 
us to determine the scaling of the diamagnetic moment 
M = Z T M R , which yields at low T: 



M(B) = — 
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plotted in Fig. [5Jd. Here Bq = %/ (2eag) is the char- 
acteristic field where scaling stops. We used the Dirac 
gas expression [l9|, [2(| for the magnetization at B (6* ) = 
Bq. From the scaling behavior of the free energy, the 
scaling of xd immeadiately follows by differentiating 
F(T, A, /x, B) twice with respect to B, while that of k fol- 
lows by differentiating twice with respect to /i. Similarly, 
one can easily derive scaling equations for mixed deriva- 
tives such as dM/d/J, or dM/dn, of interest since they are 
measurable for two dimensional electron systems [25j . 

Turning to the heat capacity, we perform the deriva- 
tive C = —Td 2 F/dT 2 and analyze the resulting scaling 
equation along the lines of our previous calculations. It 
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follows for the Dirac liquid: 
C(T) 



18T 2 C(3) 



nv 2 h 2 (1 



(12) 



The leading perturbative correction C (T) ~ T 2 — 
2AT 2 In (Tq/T) fully agrees with the recent result of 



Vafek [llj] who also argued that the dominant low-T 
behavior should be C ~ T 2 /ln 2 ^-, as follows from 
Eq. (fT2")l . At finite density, we obtain for the linear heat 
capacity coefficient of the electron or hole Fermi liquid 

7 (n) oc y/\n\/ (l + f log ^ . Using Eq. fljj]) and the 
fact that the spin contribution of the magnetic suscepti- 
bility behaves as Xs oc c (T) /T implies that the magnetic 
response at low T is dominated by orbital effects. 

Finally, we analyze the density correlation function 
Xc (<z) = (5n q Sn^ q ) where Sn q = n q — (n q ). It deter- 
mines the dielectric function e (q, u>) = 1 + ^prX (<1> w )> 
measurable in optical or electron energy loss scattering 
experiments. Due to the vanishing vertex corrections 
between fermions and collective charge fluctuations we 
find for the density n(k, u>) — Z™ n(bk,Z^ lj), yielding 
Xc = Z^b^XcR- If we use the random phase approx- 
imation result for Xc,r at b* |q| = A or Z^ 1 {b*)uo = D 7 
we obtain Xc (q, = 2ne 2 |q| /e — II (q, , but with 
the polarization (for the Dirac liquid, see also Ref [141 ]) 



n(q,w) = - 



Nq 2 



16J{-iuj + 0+) 2 +v 2 q 2 (1 + Alogx- 1 )' 



The argument of the logarithm is x 



hv(l + \logD/uj) 



and x = uj/D otherwise. 



q/A for q > 
Despite these 



nontrivial renormalizations of the charge response, it fol- 
lows once more that with a (uj) — lim^o ^jiXc (<Z, and 
the microwave conductivity of clean graphene is unaf- 
fected by interactions, as shown above. 

Our results are caused by the interaction-enhanced ve- 
locity of graphene and can be rationalized by the simple 
substitution v — > v* = v (l + Alogx -1 ) in correspond- 
ing Dirac gas expressions, where x is proportional to the 
dominating scaling variable, i.e. T, y/n, /i, y/B, q, or to. 
This insight gives an immediate physical explanation for 
the enhanced diamagnetism or the suppressed compress- 
ibility. The former is caused by the increased velocity 
of shielding currents while the latter follows from an in- 
crease in the density of states oc v~ 2 . While the simple 
rule v — > v* is trivial to implement, we stress that it is 
often the result of subtle cancellations of an infinite set of 
divergent diagrams in the perturbation theory. To deter- 
mine which scaling variable dominates for a given set of 
parameters we use standard crossover arguments. Con- 
sider for example temperature and density. For low den- 
sity, T (b) reaches its maximally allowed value T before 
n (b) reaches no and x — T/Tq, the Dirac liquid regime 
of Fig. [TJ The opposite happens at higher density, giv- 
ing x = \/\n\/riQ in the electron and hole Fermi liquid 



regimes of Fig. [TJ The crossover from one regime to the 
other takes place when n(b*) — no and T (b*) — T si- 
multaneously, yielding the crossover temperature T* (n) 



(13) 



which is shown in Fig. [TJ Similar results can be obtained 
for any pair of scaling variables. 

In summary, by exploiting the proximity to its QCP, 
we derived explicit expressions for the temperature, den- 
sity and magnetic field variation of numerous observable 
properties of graphene. This allows a direct comparison 
with experiments to reveal the role of electron-electron 
correlations in this interacting relativistic quantum liq- 
uid. 
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